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We have studied the scalar field perturbations on six-dimensional ultra-spinning black holes. We 
have numerically calculated the quasinormal modes of rotating black holes. Our results suggest that 
such perturbations are stable. 
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l/-^ , I. INTRODUCTION 

O ; 

. The production of higher-dimensional black holes in collider 0, I3| important prediction arising from 

' the large extra dimension scenario 4, 5] or the warped compactification (RSI) scenario 0, where the fundamental 
Ch (higher-dimensional) scale of gravitation can be set to the order of TeV. Such phenomena would strongly suggest the 
presence of the extra dimensions, and we would also obtain the experimental method to seek for the quantum gravity. 
Hence, it is important to investigate the fundamental properties of higher-dimensional black holes such as stability, 
which is the subject of the present paper. 

There is a variety of the higher-dimensional black objects with nontrivial topology unlike those in four-dimensional 
^ I space-times where each component of stationary black hole is a topological two-sphere. In fact, higher-dimensional 
vacuum Einstein equation admits simple solutions called black branes which are direct product of the four-dimensional 
C Schwarzschild black hole and the flat space. We are usually interested in the case where the flat dimensions in the 

black brane space-time are compactified into say a torus. 
, Gregory and Laflamme have found that the black branes are unstable under linear gravitational perturbations 
I ' when the compactification scales are large compared with the Schwarzschild radius. On the other hand, Ishibashi and 
I Kodama Q have shown that the higher-dimensional Schwarzschild black holes in the spherically symmetric space- 
times are gravitationally stable as in the case with the four space-time dimensions. An accurate way to compute 
the quasinormal frequencies of higher-dimensional Schwarzschild black holes has been obtained via 6th order WKB 
approximation 

It is also important to consider the stability of the rotating black holes, since black holes produced in the collider 
generically have angular momenta and the production cross section is larger for more rapidly rotating black holes |1C| In 
contrast to the Kerr-bound on four-dimensional rotating black holes, there are black hole solutions with an arbitrary 
large angular momentum for a fixed mass in higher-dimensional general relativity. Recently, the existence of the 
^ effective Kerr-bound on such rapidly rotating black holes in higher dimension has been conjectured by Emparan and 
■ Myers • They showed that the geometry of the event horizon of such rapidly rotating black holes in six or higher 
dimension behaves like the black membranes and argued that such rapidly rotating black holes therefore become 
unstable. 

Although we should analyze the gravitational perturbations of rotating black holes to confirm such instability, any 
formalism where the perturbative equations are separable is yet unknown. Here we just consider the perturbation 
of the massless scalar field propagating on the rotating black hole space-times. Even if we restrict ourselves to field 
perturbation on fixed background space-time, there exists an instability of massless field perturbation on rotating 
black branes and strings discovered by Cardoso and Lemos ^3 • The equation of motion for the massless scalar field 
perturbation on the rotating black holes is known to be separable in any dimension Il3l 1 141 . and authors examined 
the stability of this system in five-dimension as the simplest case in the previous work |l4| . The present paper is the 
complementary work on the renewed interest in the instability (or stability) of the rapidly rotating black holes in six 
dimensions. 
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II. MYERS-PERRY BLACK HOLE AND ITS MEMBRANE LIMIT 

The (4 + n)-dimensional Myers-Perry metric with only one nonzero angular momentum is given by 

A - sin^ 1} , , 2a{r^ + - A) . , „ , , 
g = ^ sm^ -ddtdip 

(r^ + - Aa^ sin^ i?,,!], ,,,, 

+ '—- sin^ Mlp^ + —dr^ + + cos^ dd^l, ( 1 ) 

where 

S = + cos^ (2) 

A = r'^ + - iir^-"" , (3) 

and dfi^j denotes the standard metric of the unit n-sphere. This metric describes a rotating black hole in asymptotically 
flat, vacuum space-time with the mass and angular momentum proportional to /i and /la, respectively. Hereafter, 
/i, a > are assumed. The event horizon is located at r = rn, such that A|,.=,.^ = 0, which is homeomorphic to 5'"+^. 
For n > 2, A always has a positive root for arbitrary a, namely regular black hole solutions exist with arbitrarily 
large angular momentum per unit mass. 

In the limit where a 00, the coordinate radius of the horizon is approximated by th ~ (/i/a^)^/("~^^. Hence, the 
horizon radius rn shrinks to zero when one take the limit where a — > 00 with fixed fj,. To avoid a vanishing horizon 
radius, Emparan and Myers introduced the black membrane limit where a ^ 00 and ^ —>■ 00 with fixed p, := ^/a?. 
In this limit, the Myers-Perry metric behaves like a black membrane metric near the pole d — Q: 

fl \ jj^2 , '^''^ 2 



9 = -[^-^)dt'+ + r'dni + d,^ + V, (4) 

where the new coordinate is defined by <; :— asin-d. In this paper, we are concerned with the dimensionless Kerr 
parameter a* := a/rjj. The limit where a, — > cxd with fixed rjj is consistent to the membrane limit. 



III. SCALAR FIELD EQUATIONS 

Here we consider the perturbation of the massless scalar field on the background Myers-Perry metric We put 
the scalar field configuration as 



e 



iLjt—im<f 



R{r)S{{})Y{n), (5) 



where Y{^1) is the hyperspherical harmonics on the n-sphere with the eigenvalue j{j + n — 1), {j = 0, 1, 2, • • •). Then 
we obtain the field equations with separated variables: 

( 4t sin ^? cos" ) + luj^a'^ cos^ ^ - csc^ d - j(j + n - I) sec^ iS + A]S = 0, (6) 

sint9cos"i? Vc^'? d-d J ^ j\j > j > \j 

and 

-X + iHr^+ay^aY _ j{j+n l)a^ - a1 i? . 0, (7) 

where \ := A — 2muja + uP'd} . 

We will solve the eigenvalue problem under the quasinormal boundary condition: 



where 



(n-l)(r|, + a2) + 2r|,^" 
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When a = 0, the eigenfunctions for the angular equation © are analytically given in terms of the hypergeometric 
functions: 

ft,„, = (sin79)l'"l(cos79)^i^(^-£,£ + j + |m| + ^,j + ^;cos2^) , (10) 

with the eigenvalues 

Ae,,n^{j + \m\ + 2£){j + \m\ + 21 + 71+1), {£ ^ 0,1,2, ■■ ■). (11) 

Then e^''"^'^ Sijm{^)Y (H.) behaves as the hyperspherical harmonics on 5"+^, which belongs to the irreducible represen- 
tation of the rotational group SO{n + 3). Especially for n — 2, the hyperspherical function on S"^ with the eigenvalue 
L{L + 3) belongs to the representation characterized by the Dynkin index [L, 0] and has (L + 1){L + 2){2L + 3)/6-fold 
degeneracy. When a 7^ 0, S'0(5) symmetry of the space-time breaks down to SO{2) x S'0(3) and the degenerate mode 
partially split into the several modes with (2j-|-l)-fold degeneracies, where j and m are constrained by L = j + \m\ + 2i, 
{i = 0,1,2,- ■■). 



IV. NUMERICAL COMPUTATION 



We use the continued fraction method to determine the resonant frequency lo and the separation constant A. We 
assume the following series expansion for S: 



5= (sint9)l"l(cos19)^■^afe(cos2 1?)^ (12) 



fe=0 

which automatically satisfies the regular boundary conditions at "i? = 0,7r/2 and converges at any time. Substituting 
this expansion into the angular equation 0, we obtain the three-term recurrence relations: 

atak + PkCLk-i + lkak-2 0, (fc = 1, 2, • ■ •), (13) 

where a_i = and the coefficients are given by 

au = -%k{2j + 2k + 7i~l), (14) 
= 2[{j + \m\ + 2k-2){j + \m\ + 2k + n-l)-A], (15) 
Ik = -alul- (16) 

Here a* :— a/rn and :— urn are the dimcnsionless quantities. The limit where ^ 00 corresponds to the 
membrane limit as mentioned above. 
We expand the radial function R as 

/ \ i(7 / , \ — (n+2)/2 — icT oo / \ k 



which automatically satisfies the quasinormal boundary conditions and converges at any time. When n = 2, substi- 
tuting Eq. (|17|l into the radial equation ((TJ, we obtain the eight-term recurrence relations 

akbk + I3kbk-i + jkbk-2 + Skbk^a + ekbk-4 + Ckbk-5 + Vkh-e + Skbk-7 = 0, (fc = 1, 2, • • •), (18) 
where 6_i = 6_2 = • ■ • = and the coefficients are given by 

ak = -{al + l){al+3fk{k + 2ia), (19) 
(3k = {al + 3){~2{5at + 22al + 33)a'^ -Ama^{al + 5)(7, 

+ {al + 1)[(5A:2 ~7k + 4)0^ + 3k^ - 3k + 6 + 2j{j + l)al + 2X]} 

+i{al + 3){[(14A: - 9)ai + (40fc - 22)al + 42fc - 21]cr + 2m,a4al + 3)(2fc - 1)}, (20) 
= {27al + 117at + 89al-129)(T^ + 4:ma4al + 3){5al + l)a + Am^al{al + 5) 

-{al + l)[{9k'^ - 26k + 2A)at + [Uk^ - A2k + 40)0^ - 3k^ - 24 + 4j(j + l)a*(2a^ -h 3) - 12A] 

-i{[(30fc - 42)a^ -I- (106fc - 152)a^ + (58A: - 114)a2 - 114fc -I- 108]cr + Ama^{al + 3)[(3fc A)al - 3{k - 1)]}, 
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(21) 

Sk = -{3al - l){5at + 6al - 15)a^ - 8ma,{2al + 3){al - 3)a - Am^aUal - 7) 

+ {al + l)[{5k^ - 25k + 36)at - 2{k - l)(3/c - 2)al - 3k^ + 3/c + 24 + + l)al{Sal + 2) - 4(0^ - 2)A] 
+i{[(18fc - 45)a^ + (6fc - 55)a^ + (-42/c + 25)a^ + 66fc - 141]cr + 4ma4(2fc - 5)a^ + (-4/c + 5)al + 6k - 12]}, 

(22) 

efe = -2(9a^ + 25)(a2 + l)2CT2-8ma,(2a^ + 3a^-3)cr-4m2a^(a2-3) 

+ (a^ + l)[(5fc2 - 20A: + 16)at + (UP - 60k + 64)a^ + S/c^ - 24A: + 40 - 4j(j + l)aj(2a2 + 1) + 4A] 
+i{[(18fc - 36)a^ + (78fc - 164)a^ + (llOA: - 252)a2 + 82fc - 204]cr + 4ma4(2fc - 4)a^ + (4fc - 9)al + 6k - 15]}, 

(23) 

(a2 + l){4(6at + 15a^ + 5)0-2 + 20ma»(a^ + l)a + 4m2a2 

-(fc - 3)(9A: - 28)at - (lO/c^ - 60/c + 86)al - [k - 2)(5fc - 19) + 2{al + l)[j(j + l)a^ + A]} 

-?;(a2 + 1){[(30A: - 93)at + (52fc - l6Q)al - 2k + 13]cr + 2ma4{6k - 19)0^ - 2A: + 7]}, (24) 

rik = (a2 + l)2{-(9a2 + 13)cr2-4ma,cr+(fc-4)[(5/c-18)a2 + A;-4]} 

+2i(a2 + l)2{[(7fc - 27)a2 + 7k~ 28]a + 2ma,(fc - 4)}, (25) 

Ok = (a2 + l)3[cr2-(fc-4)(fc-5)-i(2/c-9)f7]. (26) 

After we repeatedly use the Gaussian elimination as written in the Appendix, the eight-term recurrence relations H18|l 
are reduced to the three-term ones 

arbk+prh-i+irh-2 = o, (fc = i,2,...). (27) 

For given the three-term recurrence relations (|13|l . the successive coefficients are obtained in two ways, the finite 
and infinite continued-fraction representations: 

Ofe _ /?fc 1 -%dk-i -73(52 -72ai 

(28) 



a-k-i OLk ak (3k-i+ /?2+ Pi 

7fe+l -ak+llk+2 -OLk+2lk+5 



Pk+1+ (ik+2+ /3fc+3 + 

where (1 / x~\-){y / z) is an abbreviation for l/[x + {y/ z)]. Similarly for Eq. H27|) . we obtain 



bk' _ 


Pk'" 


1 ~irc^7u -ir<' 




bk'-i 












"fc' + l7fc'+2 "fe'+2 'fe'+3 








B'"" + 3""' + 

A''fe'+2 ' '''fe'+3 ' 





(29) 

For fixed k and fc', we obtain nonlinear algebraic equations (|28|l and (|29|) for two unknown complex numbers u:^ and 
A. We use a standard nonlinear root search algorithm provided by the MINPACK subroutine package to evaluate 
the numerical solutions of these equations. The equations for different sets of (/c, fc') are used as numerical check. 

When the perturbation become unstable under the membrane limit, we expect that the behavior of such instability 
resembles one of the Gregory-Laflamme instability. Since they found that the s wave is unstable, we consider the 
j = modes. The critical value of the wavelength of the Gregory-Laflamme instability is of the order of th. In the 
coordinate <j along the rotation plane, the perturbation has the wave length ^ ai? ~ a^rnl^ for £-th mode. Hence, 
we expect that the possible instabilities occur for a* > t. 

Typical examples of quasinormal modes of the six-dimensional Myers-Perry black holes are plotted in Figs.^El 
Here the real and imaginary parts of the resonant frequencies as functions of the parameter a*/-\/l + for various 
value of m and £ are plotted. Although the Kerr parameter a, is taken to have a sufficiently large value, each mode 
does not seem to have negative imaginary part in any case. 



V. SUMMARY AND DISCUSSION 



We investigated the massless scalar field equation in the six-dimensional Myers-Perry black hole background. Under 
the quasinormal boundary condition, we have searched for the resonant modes. 
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In the non-rotating case (a* = 0), the quasinormal modes with same L(= j+ \m\ + 2l) are degenerate. Continuously 
varying the rotational parameter a*, these modes split into several modes characterized by {j,m,£). The sequences 
of quasinormal modes we obtained have positive imaginary parts, thus they do not show any evidence for instability. 

Recently, Cardoso, Siopsis and Yoshida investigated the stability of the same system. Several quasinormal 
modes within our present work are also calculated by them in different way to expand the radial function and their 
results are consistent to our present results. They also conclude that the stability of scalar field perturbations on 
Myers-Perry black holes are suggested. 
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APPENDIX A: GAUSSIAN ELIMINATIONS 



The eight-term recurrence relations (|18|) are reduced to the seven-term ones 

a'kh + /3fe&fe-i + 7fc&fc-2 + + e'kh-i + Ckh-5 + v'kh-G = 0, (fc = 1, 2, • • •), (Al) 

via the Gaussian elimination: 

a[ = ai, Pi, 

a'2 = £12, = P2, 72 = 72, 

a3="3, 03 ^Ps, 73 = 73, 

a4 = a4, /?4 = P4, 74 = 74, 6'^ = 64, £4 = £4, 

"5 ="5, ^5^05, 75 = 75, S'r, = Sr-„ 4=^5, Cs = C5, 

"6 ="6, P'e^Pe, 76 = 76, (^e = (^g, 4 = <^6, Ce = C6, v'e = Ve, 
a'k = (^k, P'k^ Pk-a'^-iOk/ri'k-i, I'k ^ Ik - P'k^iOkh'k-i^ ^^I- = 4 - 7fe-i^fe/?7fe_i, 
e'k = ek-5'^,_iekh'k-i. Ck^Ck-e',_idk/Vk-i, v'k ^ Vk - CLiOk/vLi, (fc = 7,8,.-.). (A2) 
The similar procedures 



"l 


= "i. 




= /3;, 
















= "2, 


0i 


= /92, 


ii 


= 7^, 










"3 


= "3, 


P'i 


= /?3, 


ii 


= 7^, 




5'3. 








= "4, 


P'l 


= /?4, 


il 


= 7i, 




X' " 

■ 04, £4 


= 4, 






= "5' 


P'i 




1 


= 7^, 


X" 

h = 


X' " 


= 4, 


C" = 




= "fe, 


P'l 


= /3^ 


- "fe- 


-Ma- 


-1, 


7^' =7^ 






4' 


= 4- 


5'U 




-1, 


Ck = Cfe 


-4. 


-iVk/ Ck- 


-1, 


(fc = e 



^5' 



-^It = s'k-ik-iVk/Ck-i: 

••), (A3) 



/// // nl/l nil 

Oil = Oil, Pi = Pi, 

I/I II qIII nil III II 

"2 = "2 , P2 = P2 , 72 = 72 , 

/// // Qlll Oil III II x"! X'l 

"3 = "3 , P3 = P3 , 73 = 73 , = 63 , 

/// // Qlll Oil III II Cl'l X'l "' " 

"4 = "4 , P4 = P4 , 74 = 74 , f>4 = 64 , £4 = £4 , 

/// // Qlll Qll III i-ll I III III II Qlll /-Il I III X'l' X" "' /-" I "' 

Oik =Oik. Pk = Pfc - "fe-iCfe/efc-1, Ik = 7fc - Pfe-^iCfe/efc-^i, (>k ^ Ok ~ Ik-iCk / <^k-i, 

4" = 4-^1^74-1, (fc = 5,6,...), (A4) 



"1 — "1 , /Ji - Pi , 

nil III Qiiii Qlll nil III 

"2 = "2 , P2 = P2 , 72 = 72 , 

//// /// Qiiii Qlll nil III X"" X'" 

"3 = "3 ' P3 = P3 ' 73 = 73 , = C)3 , 

//// III Qlll I Qlll nil III IX"" "" "' q"" "' IX"" 

Oik =Oik, Pk Pk - "fc-iffc /'5/c-i, 7fe =7/c ~ Pk-i('k /Ok-i, 

C = C-7ri47^ri, (fc = 4,5,...), (A5) 
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FIG. 1: The quasinormal modes for {j,m,£) — (0,0,0). 



inn 

inn 
^2 
mil 



Pi — Pi : 



12 ' 



///// r//// / ///// 



Titll I lllll III. 
'>k /lk-1: Ik 



7^' 



fQiiiii X"" I mil n OA \ 

Pk-i^k /7fe-i, (fc = 3,4,---), 



(A6) 



leads to the three-term relations 



oTbk + Prbk-i + 7r&.-2 - 0, (fc = 1, 2, . . .). 



(A7) 
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FIG. 2: The quasinormal modes for {j,m,£) = (0, ±1,0). 



FIG. 3: The quasinormal modes for {j, m, £) = (0, ±2, 0) and (0, 0, 1). 




FIG. 4: The quasinormal modes for {j,m,£) = (0, ±3,0) and (0, ±1, 1). 




I 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
a/sqrt(1+a^) 



FIG. 5: The quasinormal modes for {j,m,£) = (0,±4,0), (0, ±2, 1) and (0,0,2). 




FIG. 6: The quasinormal modes for {j, m, t) = (0, ±5, 0), (0, ±3, 1) and (0, ±1, 2). 



